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Abstract 

We study the spectrum properties for a recently constructed fixed 
point lattice Dirac operator. We also consider the problem of the 
extraction of the fermion condensate, both by direct computation, 
and through the Banks-Casher formula by analyzing the density of 
eigenvalues of a redefined antihermitean lattice Dirac operator. 
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1 Introduction 



One of the most annoying drawbacks of the lattice discretization of a fermionic 
theory is the impossibility] [1 to avoid the explicit breaking of the chiral sym- 
metry. In the case of a simple discretization, the Wilson action, the breaking 
of the symmetry is so bad that no trace of the chiral properties of the con- 
tinuum theory is kept in the lattice theory. This is of particular relevance 
for attacking QCD by lattice techniques. At the classical level, the chiral 
properties of the zero-modes of the Dirac operator are lost, and the Atiyah- 
Singer theorem of the continuum theory [2j has no strict correspondence on 
the lattice. At the quantum level, the explicit breaking of the symmetry 
induces an additive renormalization of the quark mass, the chiral limit be- 
ing attained through a fine tuning of the bare parameters; moreover, chiral 
currents undergo a finite renormalization, mixing among operators with dif- 
ferent (nominal) chirality occurs, and no order parameter for the spontaneous 
breaking of the chiral symmetry can be defined in a natural way. 'Accidental' 
zero-modes (related to the so-called 'exceptional configurations') are causing 
troubles in Monte Carlo simulations. 

In an early paper [3| of lattice quantum field theory, Ginsparg and Wilson 
already gave the key to make definite what is meant by chiral limit in the 
framework of a lattice theory breaking explicitly chiral symmetry, providing 
a general condition - we will refer to it as to the Ginsparg- Wilson Condition 
(GWC) - for the fermion matrix of the lattice theory. They showed in par- 
ticular that any action fulfilling the GWC reproduces the correct triangular 
anomaly. The condition was found by requiring that the lattice action is ob- 
tained as the fixed point (FP) action of a block-spin transformation (BST), 
in the universality class of a continuum chiral-symmetric action. In this way, 
the breaking of the symmetry is introduced in the lattice action by the BST 
itself and not 'by hand'. Recently it has been shown [3] that the GWC can 
be solved in a rather independent approach, namely in the so-called overlap 
formalism, and the corresponding symmetry has been identified jH]. 

Ginsparg and Wilson's observation had no practical applications until 
the technology for the construction of FP actions of lattice gauge theory 
was improved following Hasenfratz and Niedermayer's ideas jHllHIHl- In a 
recent series of papers it has been shown that the GWC is a sufficient con- 

1 With the precondition that very general properties of the continuum theory, like lo- 
cality, are preserved. 
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dition for the restoration of the main features of the continuum (symmetric) 
theory; in particular, at the classical level, the Atiyah-Singer theorem finds 
correspondence on the lattice [HI UOj ; excluding the possibility of 'acciden- 
tal' zero-modes; at the quantum level, no fine tuning, mixing and current 
renormalization occur, and a natural definition for an order parameter of the 
spontaneous breaking of the chiral symmetry is possible [TT] . 

Monte Carlo calculations require actions with finite number of couplings, 
while the FP action is extended over all distances (even with exponential 
damping of the couplings). It is therefore clear that only a parametrized 
form of the FP action is of practical relevance. The point to be checked is 
to what extent such an approximation of the FP action is able to reproduce 
the nice properties of the 'ideal' FP action. The construction of the FP 
action of QCD is - because of technical problems - yet still far away (for 
some pioneering attempts, see ^^), and meanwhile we try to practice with 
toy models. 

Recently, a parameterization of the FP action of the Schwinger model 
for the non-overlapping BST (noBST) (we refer to this parameterization as 
to pApp subsequently, while to the 'ideal' FP action simply as to App) was 
found ^3] in a Monte Carlo approach. The fermionic action is parametrized 
in terms of bilinear expressions in the fermion fields, connected by paths made 
up of compact gauge link variables, whereas the gauge action was the non- 
compact simple fixed point action. Here we give a brief account of a study 
of the spectral properties of this Dirac operator, concerning the numerical 
verification of some of the chiral properties of the FP action in the case of the 
pApp. A more detailed report will be given elsewhere. We focus our attention 
on the spectral properties of the FP Dirac operator, verifying the lattice 
Atiyah-Singer theorem. A first verification of this theorem was accomplished 
in a different approach in ^3]. Moreover, we address the problem of the 
extraction of the fermion condensate, both from the direct computation and 
through the Banks- Casher formula [To] . 

2 Chiral Properties of the FP Action 

The GWC for the fermion matrix h xx i for massless fermions reads 

\{h x , x ,^} = (hj 5 Rh) xxl , (1) 
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where R x , x ' is a local matrix in coordinate space, i.e. whose matrix elements 
vanish exponentially with the distance. Any FP action of a BST satisfies the 
GWC: R XtX > is fixed by the average matrix of the fermion BST, and it can be 
in general assumed to be trivial in Dirac space. 

In the case of the non-overlapping BST (noBST) considered in [13] . one 
has R XjX > = ^S XjX r, and the GWC assumes the elegant form 

h + h) = h ] h = hh) . (2) 

where also the hermiticity property h> = 7 5 /ry 5 (fulfilled by the FP action in 
this study as well as by other actions) has been taken into account. 

2.1 The spectrum of the Dirac operator 

Eq. (j2J) implies non trivial properties of the spectrum of the Dirac operator, 
namely the fermion matrix h: 

i. [h, ht] = 0, i.e. h is a normal operator; as a consequence, its eigenvec- 
tors form a complete orthonormal set. 

ii. The spectrum lies on a unitary circle in the complex plane centered at 
(1,0). 

iii. The property (i), together with the hermiticity property of the fermion 
matrix, implies (denoting with v\ an eigenvector of h with eigenvalue 
A): 

l 5 v x = v x * , if A^A* 

7 5 w A = ±v x , if A = A* e R . (3) 

So, just as in the continuum the eigenvectors of complex-conjugated 
eigenvalues form chirality doublets; moreover all real-modes have defi- 
nite chirality. (For a general FP action this property holds only for the 
zero-modes.) 

It was shown in that (0) ensures the existence of a lattice version of 
the Atiyah-Singer theorem, which can be stated in the form: 

Qfp = -J2( y o^ 5v o) > ( 4 ) 

{vo} 
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where Qfp is the fixed point topological charge (THl Ej of the background 
gauge configuration. Qfp provides a lattice definition of the topological 
charge of a gauge configuration which depends on the RG blocking procedure 
for the gauge sector. 

2.2 The fermion condensate 

In a subtraction procedure for the fermion condensate - inspired by ((H) 
- was proposed: 

(^)sub = -l/y(tr(/i- 1 - J R)) gauge , (5) 

where V is the (finite) space-time volume. Because of (JIJ the quantity in 
brackets of the r.h.s of the above equation vanishes, except when a zero- 
mode of h occurs, in which case a regulator-mass fi must be introduced, 
h — > = h + fi. For the number of flavors n/ > 1, the contribution of 
the zero-modes to the gluon average vanishes when // — > because of the 
damping effect of the fermion determinant, and we conclude that the finite- 
volume subtracted fermion condensate vanishes in the chiral limit j^]. The 
case nf = 1 is peculiar since (V^sub ^ even in a finite volume. The 
configurations responsible in this case for the non-zero fermion condensate 
in a finite volume are those from the \Q\ = 1 sector; indeed, if h has just one 
zero-mode - which is possible because of the (lattice) Atiyah-Singer theorem 
only for \Q\ = 1 - the quantity tr h~ l (det h) has a non-zero limit when \i — > 0. 
The effect of the subtraction is in general (i.e. for any nf) just to remove the 
spurious contribution of the Q = sector introduced by the explicit breaking 
of the chiral symmetry of the FP action. 

The situation in the infinite volume, i.e. when the fermion condensate is 
obtained through the sequence of limits lim^ limy_+oo, is somehow different: 
in this case the role of the zero-modes is irrelevant^, the quasi-zero-modes 
being responsible for the non-zero fermion condensate as the Banks-Casher 
formula shows. 

We observe that (JSJ) may be rewritten in the form 

(^)sub = -l/^(tr(/i- 1 )) gauge , h = hil-Rh)- 1 . (6) 
The redefined fermion matrix h has the non trivial properties: 

h ] = - h , {h,l 5 } = , (7) 

2 We are grateful to P. Hasenfratz for having driven our attention on this point. 
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from which it follows that h has a purely imaginary spectrum, 

h v = -i \v , A E R; (8) 

of course, the zero-modes of h coincide with those of h. 

The replacement A — > A allows in a sense to reconstruct the 'chirally 
invariant information' of the spectrum of h. One can show that the spectral 
density of h, p(A) complies (in the limit: lim a ^ hniy^oo) with the the Banks- 
Casher formula for the subtracted fermion condensate: 

(^)sub = 

In the case of the noBST, one has 

A = A ^1 

and the spectrum of h is obtained by mapping the spectrum of h (which, we 
recall, lies on a unitary circle in the complex plane) onto the imaginary axis 
by the stereographic projection. 



-vrp(O) 



(9) 



(10) 



3 Numerical Results: the Schwinger Model 

The parametrized FP action pApp of the Schwinger model in was ob- 
tained in the non-compact formulation for the pure-gauge sector. A first 
difficulty arises in this respect, since it is not possible to give a lattice defini- 
tion of the topological charge in terms of the non-compact gauge variables. 
The most natural definition, i.e. 1/27T F 12 (x) vanishes identically on the 
torus. A definition in terms of compact variables is instead available, the 
so-called 'geometrical' charge: 

Qg = ^ Vlmln(f/ 12 (x)). (11) 

x 

Measuring the topological charge after the compactification of the gauge 
variables according to (JTTj) reveals however an unnatural suppression of non- 
zero values even at moderate /?'s, due to the non-compact nature of the gauge 
action. 
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Ideally, one should work with the FP action of the compact theory, the 
corresponding FP topological charge operator coincidinjjl in this case with 
the expression (fTTj) [T4*] . Since a parameterization of a compact FP action is 
not yet available, here, as an approximation, we take the parametrized FP 
fermion matrix of ^3] (which is compact by construction) and simply replace 
the original action of the pure-gauge sector with the compact Wilson action. 
Of course this approximation is expected to introduce additional deviations 
from the behaviour of the ideal FP action App when topological fluctuations 
are important in the statistical ensemble. 

The fermion part of pAp p has the form 

3 

^ h p (U) V> = ^ J>i(/) $(x) <n U(x, f) ij(x + 5f) . (12) 

i=0 xj 

Here h p (U) is the parametrized lattice Dirac operator, / denotes a closed 
loop through i or a path from the lattice site x to x + 5f (distance vector 
Sf) and U (x, f) is the parallel transporter along this path. The <7j-matrices 
denote the Pauli matrices for % = 1,2,3 and the unit matrix for i = 0. We 
make the identification: 7 = a±, 7 1 = 02, 7 5 = cr 3 . The action obeys the 
usual symmetries as discussed in |T3] ; altogether it has 429 terms per site. 

We constructed gauge field configurations according to (a) the compact 
Wilson action S\y and (b) the non-compact action S^c — (/5/2) J2 x ^i2( x )- 
In the case (a) the configurations were generated with a Metropolis Monte 
Carlo update, separating configurations by a number of updates of twice 
the size of the integrated autocorrelation length for the (geometric) topo- 
logical charge, which scales with (3 exponentially (e.g. approximately tq ~ 
exp (1.67/5 — 3) for (3 > 2 on 16 2 lattices; for comparable observations see 
18J). The measured distributions for the topological charge agree with the 
observations in [T9*] . In situation (b) the configurations were generated com- 
pletely independently according to the Gaussian measure (respecting the 
gauge d.f.). 

For each of those gauge field configurations we analyzed the Dirac opera- 
tor eigenvalues (and sometimes eigenvectors) and the determinant. So most 
of the results presented may be considered quenched results, although some 
of the observables have been weighted with the determinant. The latter men- 
tioned numbers are for the full, dynamical system, although maybe plagued 

3 Except for a set of measure zero, when Ui2(x) = — 1 for some x, in which case the 
prescription Ijlljl is not well-defined. 
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Figure 1: Eigenvalues of the lattice Dirac operator (jl2J) at the values of 
(3 = 2,4,6 (from left to right), on 16 2 -lattices and sampled over 25 gauge 
configurations each, according to the compact (upper row) or non-compact 
(lower row) gauge action. 

by the expectation value of the determinant in the denominator - as we will 
discuss below. 

3.1 The spectrum 

Fig. [T] shows the spectrum of the studied lattice Dirac operator (|T2^) col- 
lectively for 25 gauge configurations (uncorrelated to the amount discussed 
earlier; we collected data for several hundred such configurations, but do not 
plot them in order to prevent the figure files becoming too large. The overall 
behaviour is well represented in the figures). At larger values of (3 we notice 
excellent agreement with the expected circular shape for App- For smaller /3 
the spectrum roughens, although the circular shape generally is kept. 

One has to keep in mind that the numerically determined parametrized 
Dirac operator (|12J) is truncated and optimized on a finite sample of (non- 
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compact) gauge field configurations. Thus it is not surprising, that the spec- 
tral shape is closer to the optimum for the configurations sampled according 
the non-compact gauge action in Fig. ^ 

Subsequently we discuss (except when explicitly stated differently) our 
results for the compact gauge action, for rif = 1 and lattice volume 16 2 . 

We note, that there are configurations with real eigenvalues. We checked 
the eigenvectors for those and confirm that these modes have definite chiral- 
ity, indeed, as expected from (jSJ). Also, we can clearly distinguish the real 
values around zero from those around 2 (right-hand part of the spectrum). In 
fact each real quasi-zero-mode has three real 'doubler'-partners at the other 
end of the spectrum. The chiralities of the partners add up to the opposite 
value of the chirality of the quasi-zero-mode(s). 

We may identify these real eigenvalues (around zero) with zero-modes 
and relate their number no with the geometrically (i.e. from the gauge field 
configuration) defined topological charge Qq. We find agreement in the fol- 
lowing sense: The ratio of the number of configurations, where these numbers 
coincide over those, where they do not, approaches unity in the limit (3 — > oo. 
In fact, the approach is faster than that observed in the situation of the Wil- 
son Dirac-operator [TH]. We find that ratio to be 0.887 , 0.999 and 1.000 at 
P = 2, 4 and 6, respectively. 

Whereas App has no real modes at negative values, here we do observe 
those, in particular towards smaller (3. This is expected due to the general 
roughening (cf. Fig. [TJ) and is equivalent to the (albeit more frequent) occur- 
rence of 'exceptional' configurations for the Wilson fermion action already 
below k c . 

In order to estimate the scaling behaviour of the deviations of the spec- 
trum from the ideal circular shape, we defined a mean deviation |A — 1| 
from the unit circle in the region close to A = 0, in an angular window of 
|arg(l— A) | < 7r/4. Fig. Elshows the behaviour of the average width (standard 
deviation) o of that distribution with regard to (3. The log-log plot demon- 
strates a behaviour of a oc 1/ f3 2Al ~ a 5 . The parametrized action pApp is 
truncated in a finite range (7x7 in our case). Heuristically this implies an 
error for the eigenvalues, considered as dimension-one gauge-invariant oper- 
ators of the gauge field, in the form of some operator of higher dimension k. 
From the observed deviation we estimate an effective value k ~ 5. Since sim- 
ilar behaviour (differing just in an overall multiplicative factor) is observed 
for both types of gauge action, we think that we may justify the observed 
deviation by the truncation. 
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Figure 2: Scaling behaviour of the distribution width a for the spectra 
observed in Fig. 1 a-c (full circles, compact gauge action) and Fig. 1 d-f 
(diamonds, non-compact gauge action) as discussed in the text. The lines 
are the result of a linear fit. 

3.2 The fermion condensate 

Spectral density: The spectral density for the projected eigenvalues A 
(see (JIHl)), p(A), has been calculated both in the quenched and unquenched 
situation. In the unquenched case, the individual configurations are weighted 
by the value of the determinant, which amounts to sampling in the presence 
of dynamical fermions. The fluctuations are in this case higher, since one 
has to divide by the sum over determinant values, and as is well known 
this involves cancellation of large terms and therefore substantial statistical 
errors. 

The results are displayed in Fig. El for a sample of 2000 configurations. 
In the quenched case we observe in the central part of the spectrum a peak 
(smoothed by the binning procedure) produced by the (quasi-) zero-modes. 
It is removed when they are discarded from the sample (cf. Fig. 0J). This 
peak is absent in the unquenched situation as expected: The zero-modes are 
suppressed due to the determinant. Apart from this, the unquenching does 
not change the main features of the distribution, only introducing higher 
fluctuations. For both, the so 'cleaned' distribution as well as for the full, 
unquenched case, we observe good agreement of p(0) with the theoretically 
expected value for the condensate at infinite volume (ipip) = —e 1 /(2iry/j37T) 
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Figure 3: Density of the eigenvalues (projected according (fTUj) from the 
close-to-circular positions to the close-to-imaginary axis positions) near zero; 
(3 = 4, volume 16 2 . The horizontal line gives the theoretical expectation from 
the Banks-Casher formula. Left: Quenched histogram and histogram without 
zero-modes (dotted). Right: The histogram weighted by the determinant, 
i.e. the unquenched distribution. 



through the Banks-Casher relation (JHJ). The theoretical value in the finite 
volume case (torus) is also available (201; for the physical volume here con- 
sidered the correction is less than 4%. 

Our conclusions from these observations are: 

• The partial unquenching (by omitting the zero-modes) gives similar 
results as the complete unquenching (by including the determinant). 

• The zero-modes produce a substantial finite-volume effect; indeed their 
elimination allows us to approximately recover the Casher-Banks rela- 
tion even in the quenched case. 



Direct computation: Alternatively we may attempt to derive {4>ip) di- 
rectly from tr(/i _1 ) according to (jSJ. 

The standard procedure to get the infinite volume fermion condensate 
is to take the succession of limits lim^o limy^oo. We argue that we have 
indeed accomplished something equivalent to this in the previous calculation 
with the spectral density, via the Banks-Casher formula: For that quantity, 
the volume effects are small (after discarding quasi zero- modes), so that the 
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Figure 4: We show —(ipip) S ub,fi determined from the same sample of config- 
urations as discussed in Fig. 01 

V — > oo limit is effectively attained, and the limit A — > replaces the limit 
/i — > 0. Note, that these two limits are complementary, since \i represents 
the real part of the eigenvalues of h, while A the imaginary part. One can 
show j2Ij that the two approaches, the direct computation and the one via 
the Banks-Casher formula, are consistent even in a finite volume, if \i ^> 
1/V(fi(>). ' 

First of all, we have checked that, as expected from the discussion of 
the previous section, the subtracted fermion condensate © is zero (within 
error bars) in the trivial topological sector. We then studied the full con- 
densate, i.e. including also (quasi-) zero-modes. Fig. |U gives — (^V) sub,^ for 
the discussed sample of 2000 configuration (/3 = 4, volume 16 2 ) determined 
for several values of the regulator mass fi. The values are of course strongly 
correlated since each point contains the same set of eigenvalues of h. The 
comparison with the (infinite volume) continuum value shows good agree- 
ment near /i = 0. In fact, it turns out, that the regulator here is not really 
important. Essentially all quasi- zero- modes are not exact zero modes and so 
the regularization of the determinant (vanishing for an exact zero mode) and 
the inverse eigenvalue (diverging) is not necessary. 

Obviously a better study of the finite volume dependence for the presented 
quantities is in order. However, already from the results presented here we 
conclude, that - at least for the rif — 1 Schwinger model - even approximate 
fixed point actions provide excellent possibilities to study chiral properties. 
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